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Recently we have predicted [Phys. Rev. Lett. May 2005] that Ti-decorated carbon nanotubes 
can absorb up to 8-wt% hydrogen at ambient conditions. Here we show that similar phenomena 
occurs in light transition-metal decorated Ceo- While Sc and Ti prefers the hexagon (H) sites with 
a binding energy of 2.1 eV, V and Cr prefers double-bond (D) sites with binding energies of 1.3 and 
0.8 eV, respectively. Heavier metals such as Mn, Fe, and Co do not bond on Ceo- Once the metals 
are absorbed on Ceo, each can bind up to four hydrogen molecules with an average binding energy 
of 0.3-0.5 eV/H2. At high metal-coverage, we show that a C60 can accommodate six D-site and 
eight H-site metals, which can reversible absorb up to 56 H2 molecules, corresponding to 7.5 wt%. 

PACS numbers: 61.46.+w,68.43.-h,84.60.Ve,81.07.De 



An efficient storage media for hydrogen is crucial for 
the advancement of hydrogen and fuel-cell technologies 1 . 
There have been a great number of reports on the 
search for new routes to engineer nanomaterials so that 

(a) they dissociate H 2 molecules into H atoms and 

(b) reversibly adsorb hydrogen molecules at ambient 
conditionsiiSii^iiSiSiLSiii. Much effort has been focused 
on the engineering of carbon-based materials such as 
nanotubesi2*ii and metal hydrides such as alanatesi^. It 
is found that while hydrogen-carbon interaction is too 
weak 10 , the metal- hydrogen interaction is too strong for 
hydrogen storage at ambient conditions. Very recently 
we have showniii a novel way to overcome this diffi- 
culty by forming artificial metal-carbide like structures 
on carbon singled- wall nanotubes (SWNT). From accu- 
rate first-principles calculations, we show that a single 
Ti-atom adsorbed on a SWNT can strongly bind up to 
four hydrogen moleculesi 3 -. At large Ti coverage we find 
that a (8,0) SWNT can store hydrogen molecules up to 
8-wt%, exceeding the minimum requirement of 6-wt% for 
practical applications. Even though these results were to- 
tally unexpected, we explained them by a simple Dewar, 
Chatt, and Duncanson (DCD) modeli^i^, where the in- 
teraction is due to donation of charge from the highest 
occupied orbital of the ligand to the metal empty states 
and a subsequent back donation from filled d-orbitals into 
the lowest unoccupied orbital of the ligand. 

Here we show that similar phenomena also occurs in 
light-transition metal decorated Ceo molecules. Below 
we first discuss several possible absorption sites for a sin- 
gle Ti atom on a Cgo molecule. We then show how a 
single Ti atom on a Cgo can bind up to four hydrogen 
molecules via Kubas interactioniii^. Multiple metal cov- 
erage cases, yielding up to 8 wt% hydrogen absorption, 
are discussed next. Finally, we briefly discuss the results 
for other transition metals from Sc to Co. 

The energy calculations were performed within the 
plane-wave implementation 16 of the generalized gradi- 
ent approximation^ to DFT. We used Vanderbilt ultra- 
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FIG. 1: A single Ti atom absorbed at hexagonal (H) (a) and 
pentagonal (P) (b) hollow sites, and double (D) and single 
(S) bond sites of a C60 molecule, respectively. 



soft pscudopotcntialai^ treating the following electronic 
states as valence: Ti: 3s 2 3p 6 3d 2 4s 2 , C: 2s 2 2p 2 and H: Is. 
The cutoff energy of 350 eV is found to be enough for to- 
tal energies to converge within 0.5 meV/atom. The calcu- 
lations are carried out in a cubic supercell of 16 A of side 
for single-metal Cgo systems and 20 A for full-coverage 
cases. We also carried out spin-polarized calculations for 
cases where the ground state of the metal-coated Ceo is 
magnetic. Structural relaxations were considered to be 
converged when forces on atoms were smaller than 0.2 
eV/A. 

Figure n shows the four possible absorption sites on a 
Ceo molecule that are considered in this study. The Ceo 
molecule has the truncated-icosahedral symmetry with 
twenty hexagons (H), twelve pentagons (P), thirty dou- 
ble C-C bonds (D) between two hexagons and sixty single 
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C 60 Ti(H) 


C 60 Ti(P) 


CeoTi(D) 


CeoTi(S) 


d(Ti-C) (A) 


2.27 


2.33 


2.10 


2.24 


d(C-C) (A) 


1.42/1.45 


1.44 


1.52 


1.48 


Q(Ti) (e) 


1.39 


1.09 


0.99 


0.84 


S(Ti) 


0.99 n 


1.40 n 


1.13 h 


1.54 h 


E s (Ti) (eV) 


2.098 


1.633 


1.837 


1.220 



TABLE I: The calculated Ti-C and C-C bond distances, the 
Mulliken charges, spins and the binding energies for a single 
Ti atom absorbed at four different sites of a C60 molecule as 
shown in Fig0 For bare Ceo the calculated double and single 
bond lengths are 1.44 A and 1.38 A, respectively. 



C-C bonds (S) which are between pentagon and hexagon 
carbon rings. We calculated the binding energy of a sin- 
gle Ti atom at these sites, which are listed in Table H] 
along with relevant structural parameters. The binding 
energy is defined as 

Eb (Ti) = E(C 6Q ) + E spol {Ti) - E spol {C 60 Ti), (1) 

where E spo i is the spin-polarized energy. Hence the posi- 
tive binding energy indicates the stability of the system. 

The results in Table [I] indicate that H-site is the most 
stable configuration for CgoTi with a binding energy of 
2.1 eV. The D-site comes next with the shortest Ti-C 
bonds among the four possible configurations. For all 
cases, we have about one electron charge transfer to the 
Ceo molecule. The S-site is the least stable absorption 
site and therefore we do not consider it any further. 

Figure 2(a) shows the energy variation from spin- 
unpolarized calculations as a single H2 molecule ap- 
proaches CeoTi(D). The energy first decreases slowly as 
the hydrogen gets closer to the C60 and Ti. However, 
as the charge overlap gets large, the H 2 molecule is at- 
tracted towards the Ti atom with a sudden decrease in 
the energy. At this point, the H2 molecule is still in- 
tact with a significantly increased H-H bond length of 
0.9 A. The second sudden decrease in energy is achieved 
by dissociating the H2 molecule into two H atoms. At 
this point, the H-H distance increases from 0.9 A to 2.94 
A. The interaction between H2 and CeoTi(D) is always 
attractive and therefore H2 is absorbed onto a Ti atom 
without any energy barrier. The final geometry is shown 
in the inset to Fig. 2(a), with relevant structural param- 
eters giveniS. In order to calculate the binding energy 
for this dissociative adsorption, we calculate total ener- 
gies of the initial CgoTi(D) and H2 states and the final 
C6oTi(D)H2 state (dashed lines in Fig. 1(a)) from spin- 
polarized calculations. We obtained the binding energy 
to be 1.16 eV (Fig. 2(a)). 

Figure 2(b) shows the energy variation as two hydro- 
gen molecules approach the Ti atom; one from each side 
of the TiH 2 group. As in the case of single adsorp- 
tion, the energy always decreases, first slowly and later 
very rapidly at which point both hydrogen molecules are 
strongly attached to the CeoTi(D)H2 system. In the final 
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FIG. 2: Energy versus reaction paths for successive dissocia- 
tive and molecular adsorption of H2 over a single CeoTi(D). 
(a) H 2 + C 60 Ti(D) -^C 60 Ti(D)H 2 . (b) 2H 2 + C 60 Ti(D)H 2 -> 
C 60 Ti(D)H 2 -2H 2 . (c) H 2 + C 6 oTi(D)H 2 -2H 2 -» C 60 Ti(D)H 2 - 
3H 2 . The zero of energy is taken as the sum of the energies 
of two reactants. 



configuration, we note that two H 2 molecules were ro- 
tated by 90° as shown in Fig. 2(b). We denote the final 
product as C6oTi(D)H 2 -2H 2 . The total energy change 
upon adsorption is about 0.65 eV (i.e. 0.325 eV/H2). 
Unlike the first adsorption, the two hydrogen molecules 
are in intact but with a rather elongated bond length of 
0.79 A. 

Figure 2(c) shows the energy evolution when a fourth 
hydrogen molecule approaches the C6oTi(D)H 2 -2H 2 sys- 
tem from the top. The energy again decreases continu- 
ously, indicating a zero-energy barrier. The final product, 
denoted as C6oTi(D)H 2 -3H 2 , is shown in the inset with 
the relevant structural parameters. The energy gained 
by the fourth adsorption, 0.433 eV/H2, is slightly larger 
than for the previous one. The H-H distance of the top 
H2 is 0.85 A. Several attempts to add a fifth hydrogen 
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FIG. 3: Molecular absorption of a single (a) and four H2 
on CeoTi(H). The binding energies and relevant structural 
parameters are also given. 



molecule at a variety of positions failed, suggesting a limit 
of 4 H2 /Ti. However in the view of a recent study^S which 
shows that it is possible to attach twelve hydrogen atoms 
to a single transition metal, it is quite possible that there 
may be other transition paths that could yield more than 
four hydrogen molecules per CgoTi. 

We next discuss the hydrogen absorption properties of 
a Ti on the H-site of a Ceo as shown in Fig. 3. For the 
first hydrogen absorption, we find that the H2 stay in- 
tact unlike the case of Ti on a D-site where the first H2 
is absorbed dissociatively as shown in Fig. 2(a). This is 
due the fact that Ti is absorbed on H-site very strongly 
and therefore there is not enough charge left on Ti to 
transfer to the a* of H2 to dissociate it 13 . The binding 
energy is about 0.58 eV and H-H bond length is 0.813 
A. Introducing additional hydrogen molecules yield more 
molecular absorption without any energy barrier up to 
4H2/Ti, which is denoted as C6oTi(H)-4H2 and shown 
in Fig. 3(b). The final optimized structure is very sym- 
metric, all hydrogen molecules benefit equally from the 
bonding with Ti atom. The average binding energy per 
H2 is about 0.465 eV, slightly smaller than the one of the 
first absorption. We have also calculated the binding en- 
ergy for the isomer, CeoTi(H)H2-3H2, where the first H2 
is bonded to Ti dissociatively. We find that this isomer 
is about 0.1 eV higher in energy than CgoTi(H)-4H 2 . 

The results for the absorption of hydrogens on the 
P-site Ti are summarized in Fig. 4. The first absorp- 
tion is found to be dissociative without activation en- 
ergy. The binding energy is about 0.615 eV, significantly 
smaller than that of Ti on D-site (see Fig. 2(a)), but the 
absorption is still dissociative unlike Ti on H-site (see 
Fig.3(a)). The two isomers, C 60 Ti(P)H 2 -3H 2 (Fig. 4(b)) 
and C 60 Ti(P)-4H 2 (Fig. 4(c)), are found to be almost 
degenerate and yield a binding energy of 0.43 eV/H 2 . 

To this point we have discussed the interaction of H 2 
with a single Ti atom bonded to a Cgo, but clearly one 
can imagine attaching more Ti to a Cgo, thereby increas- 
ing the hydrogen storage capacity. In order to show the 
feasibility of this approach, we consider several cases. 
Fig. 4(d) shows the full coverage case where all the P- 
sites of a Ceo molecule is absorbed by Ti atoms, each 
of which binds four H 2 molecules. The calculated aver- 



(a) C 60 Ti(P)H 2 (b) C 60 Ti(P)H 2 -3H 2 




(c) C S0 Ti(P>4H 2 (d) C^TKPMH,],, (6-5 wt% Hj) 




E„= 1.718 eV/4Hj E„= 0.4921 eV/Hj 



FIG. 4: The optimized structure of one (a) and two isomers 
(b-c) of four hydrogen molecules absorbed on a single Ti atom 
on P-site of a Ceo- (d) Full coverage case where all twelve 
pentagon sites are decorated by Ti, each of which absorbs 
four H2 molecules. 



age binding energy for each hydrogen molecule is 0.492 
CV/H2. This is slightly higher than the single-coverage 
case as shown in Fig. 4(c). This is due to the more sym- 
metric configuration of the full coverage case, where Ceo 
does not have to distort locally as in single-coverage case. 
In fact, we see the same effect in the binding energy of 
Ti atom on P-site, which is 1.633 eV (see Table HJ, for a 
single Ti atom. However it becomes 2.115 eV/Ti when 
all the P-sites are absorbed by Ti. These results are 
quite promising, suggesting that the full coverage system 
is more stable than the single-atom coverage. 

Figure 5(a) shows another case where six Ti atoms ab- 
sorbed on the D-sites (those sites along the two-fold axis 
of the Ceo molecule) . As in the case of full pentagonal- 
coverage, the average binding energy per H2 is slightly in- 
creased to 0.592 eV from the single Ti value of 0.559 eV. 
Interestingly, one can further add eight more Ti-atoms at 
the hexagonal faces of the Ceo molecules (those along the 
[111] directions), yielding total 14 Ti atoms per C§$. The 
full hydrogen absorbed case of this coverage is shown in 
Fig. 5(b). The average binding energy is 0.522 eV/H2, 
in excellent agreement with the 0.505 eV/H2 based on 
the single Ti-atom cases. This indicates that the 14 Ti 
atoms and 56 hydrogen molecules shown in Fig. 4(b) are 
not too close to each other and therefore suggest that the 
system has the capacity to have many Ti and hydrogen. 
In fact, the configuration shown in Fig. 5(b), which have 
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Properties 


C 60 Sc(H) 


C 60 Ti(H) 


C 60 V(D) 


C 60 Cr(D) 


C 60 Mn(D) 


C 60 Fe(D) 


C 60 Co(D) 


d(TM-C) (A) 


2.29 


2.27 


2.18 


2.09 


2.30 


2.28 


2.26 


d(C-C) (A) 


1.42/1.46 


1.42/1.45 


1.46 


1.50 


1.45 


1.45 


1.45 


Q(TM) (e) 


1.55 


1.39 


0.78 


0.91 


0.97 


0.86 


0.76 


S(TM) 


0.24 n 


0.99 h 


2.0 h 


2.35 n 


2.99 h 


2.53 % 


1.54 h 


E S (TM) (eV) 


2.127 


2.098 


1.308 


0.760 


-0.017 


-0.130 


-0.503 


E S (H 2 ) (eV) 


0.300 


0.454 


0.497 


0.239 


-0.092 







TABLE II: The calculated TM-C and C-C bond distances, the Mulliken charges, spins and the binding energies for a single 
TM atom (TM=Sc, Ti, V, etc) absorbed on a Ceo molecule. 



(a) C 60 [Ti(B)H 2 -3IL,] 6 (l>) C 60 [Ti(B>H 2 -3H 2 ] 6 [Ti(H)-4H 2 ] 8 




4.6 (3.4)wt%H 3 • 

E B = 2.368 eWH 2 -3H 2 (2 235 eV) 



7.5 (6.7)wt%H 3 o 

^=2.089 eV/4H 2 (2.021 eV) 



FIG. 5: Two high-density coverage on Ti-coated Ceo, cor- 
responding to 4.3 wt% and 7.5 wt% hydrogen storage when 
all hydrogens can be released and 3.4 wt% and 6.7% storage 
when only molecular ly absorbed hydrogens are released. 



the chemical formula, CeoTii4 H56, stores approximately 
7.5-wt% hydrogen. 

It is important to know if the results reported above for 
the CgoTi system hold for other transition metals. There- 
fore we have also studied the transition metals from Sc 
to Co. Table |n] summarizes our results for the binding 
energies and relevant structural parameters. Briefly we 
find that while Sc and Ti prefers the H-site, V and Cr 
prefers the D-sites. The binding energies monotonically 
decreases as we move from left to the right of the peri- 
odic table. In fact for heavy transition metals such as 
Mn, Fe, and Co, we obtained negative binding energy, 



indicating that the system is not stable against Ceo and 
isolated metal atom formation. We also confirmed that 
similar hydrogen absorption happens for the other transi- 
tion metals. The last line in Tablc[n]mdicates the average 
binding energy per H2 in C60TMH2-3H2 configurations 
where TM=Sc, V, etc. It increases from 0.3 eV/H2 for 
Sc to 0.5 eV/H 2 for V and then decreases to 0.239 eV/H 2 
for Cr and then the system becomes unstable for heavier 
transition metals. 

In conclusion, using the state-of-the-art first-principles 
calculations we show that light-transition metal- 
decorated Cgo molecules exhibit remarkable hydrogen 
storage properties. Metals bonded on the D-site disso- 
ciate the first H2 molecule without any activation barrier 
and then reversible and molecularly bond three more hy- 
drogens. However metals absorbed at the H-sites binds 
four hydrogen molecules without dissociating it. Com- 
bining these two binding sites we shows that a single 
Ceo molecule is able to bind up to 56 H2 molecules, 
corresponding to 7.5 wt% storage capacity. These re- 
sults combined with our previous work in Ti-decorated 
SWNT - suggest a new direction to high-capacity hy- 
drogen storage materials by decorating light-transition 
metals on nanostructured materials. The nature of hy- 
drogen bonding is explained by DCD modeli^i^ and has 
the right strength for room temperature reversible hy- 
drogen storage. 

This work was partially supported by DOE under 
Grant No. DEFC36-04-GO14280. 
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